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ABSTRACT 


Laminar flow heat exchangers which cool oil in non-inter- 
connected parallel passages can experience non-uniform flows and 
a reduction in the effective heat exchanger coefficient in a range 
of Reynolds number which varies with tube length and diameter, 
tube wall temperature and fluid inlet temperature. The method of 
predicting the reduction in effective heat transfer coefficient 
and the range of Reynolds number over which these instabilities 
exist is presented for a particular oil, Mobil aviation oil 120. 
Included, also, is the prediction of the effect of radial viscosity 
variation on the constant property magnitudes of friction and heat 
transfer coefficient. 
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INTRODUCTION 


Heat exchangers that operate in the laminar flow regime with non- 
interconnected parallel passages can depart from predicted performance as 
a result on non-uniform flow distributions. Variations in the flow dis- 
tribution among the tubes can be caused by the effect of the fluid vis- 
cosity/temperature relation, natural convection and non-uniform passage 
sizes (from design and manufacturing tolerances). If the exchanger is 
oriented vertically with the colder side up, the effect of natural con- 
vection can be eliminated. Non-uniform passage size can be minimized 
through quality control. Therefore, the effect of the viscosity/tempera- 
ture relation on non-uniform flow will be analyzed in this paper. 

The viscosity/temperature relation for the Mobil aviation oil 
120 used in this analysis is presented in Fig. 1. Viscosity is a strong 
function of temperature, decreasing as temperature increases. In this 
analysis it is assumed that other fluid properties such as density and 
thermal conductivity are uniform at 90 = 55.0 mone (881.0 ey hak 
k = 0.08 BTU/hr-ft-F (0.138 W/m*C) and cp = 0.5 BTU/Ibm-F (2.09 x 10° 
J/kg-C ). This is a reasonable assumption for oils. Also, it is assumed 
that pressure drop entrance and exit effects can be neglected. The 
analysis was performed for flow in tubes. Uniform wall temperature was 
ussed assuming the heat transfer coefficient was very much greater on 
the outside than on the inside of the tube. 

When cooling liquids or heating gases, non-uniform flow can 


exist in the flow passages. For the same overall pressure drop, two 





discrete flow rates can exist in the same heat exchanger. 
This problem has been discussed by Mueller [1], [2] and [3] 


and Rohsenow [4]. 


EFFECT OF RADIAL VISCOSITY VARIATION 


The fully developed fluid and temperature distributions in a 
tube are non-parabolic due to the viscosity/temperature relation. This 
causes a lower velocity of the colder fluid near the tube wall and a 
higher velocity of the hotter fluid near the centerline compared with 
the constant property case (ie. uniform viscosity) with the same flow 
rate and bulk temperature. 

For steady state, fully developed flow and viscosity as a 


function of temperature, the momentum equation reduces to: 


dp ia ou 
ro ene (ru ae, (1) 


Neglecting the effects of axial conduction and assuming constant heat 
flux at the tube wall, the temperature distribution for fully developed 


conditions is: 


oT d 
CE ge = 2. [ Cp u Trdr (2) 
An expression is needed to account for the change in friction 
factor and Nusselt number for the variable properties case (ie. variable 


viscosity) compared with the constant properties case (ie. uniform 





viscosity). Assume that the effect can be expressed as the ratio of the 


wall to bulk viscosities raised to an exponent as follows: 


f 
(& 


’s) ° 





p 
n 
Nu _ Mo 
Nu ee os) 
C Mb 
p 
f. and Nu. are based on constant properties at the bulk 
p P 


temperature. 

To find exponents m and n, Eqs. (1) and (2) were solved simul- 
taneously using an iterative procedure. An average fluid velocity, wall. 
temperature and parabolic temperature distribution were assumed initially. 
The fluid viscosity versus temperature relation of Fig. 1 was closely 
approximated by a piecewise curve fit using a sixth degree polynomial. 
Based on the assumed temperature distribution, the viscosity u(r) was 
determined. Equation (1) was solved for the fluid velocity based on the 
viscosity u(r). The dp/dx term was adjusted to obtain the desired 
average fluid velocity. Next, the temperature distribution from Eq. (2) 
was determined using the calculated fluid velocity. Then, the friction 
factor and Nusselt number were calculated. The viscosity u(r) was de- 
termined for the revised temperature distribution and the procedure was 
repeated until the calculated friction factors and Nusselt numbers agreed 
within less than 1%. 


The results of this analysis showed that the exponents m_ and 





n varied as a function of wall and bulk temperature as shown in 


Pies. 


Zand 3: 


fluid velocity. 


SINGLE TUBE PRESSURE DROP 


The exponents were invariant to changes in the average 


Having quantified the radial effects of the viscosity/tempera- 


ture relation on friction factor and Nusselt number using the ratio of 


the wall to bulk viscosities and the exponents 


is performed on a heat exchanger tube. 


m and n, an analysis 


An energy balance at the tube 


wall for uniform wall temperature yields: 


or 


we dT = ht D(T - T )dx 
p O 





Preliminary calculations showed 


(5) 


(6) 


that the fluid flow was fully 


developed a short distance from the tube entrance and entrance effects 


could be neglected. 


However, entrance effects on the temperature dis- 


tribution had to be taken into account because of the slowly developing 


temperature distribution. 


The following equation was suggested by 


Kays [5] to account for variation in the Nusselt number in the entrance 


region of tubes for uniform wall temperature: 


4 “Cp 
0.0668 tx k 


Nu = 3.66 + ——-+ “*_*___ ), 


ie 


~<P | 
1 +0.04 on Pa 


=} 


x 


(7) 
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Therefore, the equation to determine the temperature drop along the heat 


exchanger is given by combining Eqs. (4), (6) and (7): 


“4 , we 
AV ee ae UdToste Tx k Tk 
O b flere 0.04 [-— =a p 


Equation (8) includes the radial effects and the developing temperature 
distribution effects. The temperature change along the tube was solved 
by iteration of values of bulk viscosity, local temperature and exponent 
n. The temperature change and bulk viscosity along the tube for various 
flow rates is shown in Figs. 4 and 5 for a 0.25" diameter tube and tube 
wall temperature of 40 F. 

It can be shown that the pressure drop for an increment dx 


along a tube is given by: 


1 G2 
~ Re D 20 0) 


a8 
af 


Equation (9) is multiplied by Eq. (3) to account for radial effects on 


the friction factor. 


m 


dx Re, Hy D 20 


In integral form Eq. (10) becomes: 





al 


Ec 
mess) m  il-m 
AP = To D4 Ww J us dx (11) 


The curves of pressure drop versus flow rate for various lengths 
of tube were calculated (Fig. 6) for D = 0.25", qT) = 40 F and T, = 300 F. 
It is observed that AP goes through a maximum, then a minimum, followed 
by a continuous rise as the flow rate is increased. From Eq. (11), 
AP ~ ee therefore AP vs. w curves for other tube diameters would 
look like those of Fig. 6 but be displaced vertically proportional 


EO pv. 


MULTI-TUBE HEAT EXCHANGER PERFORMANCE 


These curves are for single tubes but may be used to predict 
the flow maldistribution and resulting effective Nusselt number for 
multi-tube heat exchangers. This exchanger performance is determined 
by how the total flow is supplied - constant flow pump cee pump), 


constant head source or variable head-flow pump (centrifugal pump). 


a) Constant Flow Pump 


In the case of a heat exchanger with five tubes (L = rime 
D = 0.25"), non-uniform flow can exist as shown in Fig. 7. For an in- 
let flow rate of 25.75 lbm/hr (average flow rate of 5.15 lbm/hr per 
tube) and inlet temperature of 300 F and wall temperature of 40 F, 
four of the tubes have a flow rate of 1.26 lbm/hr with an exit tempe- 


rature of 42.7 Ff and one tube has a higher flow rate of 20.7 lbm/hr 
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with an exit temperature of 220.0 F. The average bulk temperature at 
the exchanger outlet is 185.2 F. Figure 8 shows a plot of temperature 
along the exchanger for one tube at higher flow and four tubes with 
lower flow. The average exit temperature (Tayg) is the summation of 
the individual tube exit temperature weighted by the percentage of 
flow through each tube. This is the temperature that would be mea- 
Sured at the heat exchanger exit after mixing the five streams. If 

it were not known that this maldistribution existed and it were as- 
Sumed that the flow was uniform in the five tubes at one-fifth of 

the total flow, the effective Nusselt number would be tetermined 


using Taveg as follows: 


N = Rege? 
Veft k 
(12) 

c W Ty - qT; 

Tape 7 = qo, ole (5 - T 
O O ave 

For this case h = 5./2 BTU/hr-ft- °F (32.49 Svea and Nu = 1.49. 

eff ene 


Also shown in Fig. 8 is the temperature along the exchanger and the exit 


temperature (T ) if the tubes, in fact, had uniform flow. For 


uniform 
; 4 2. ig : 
this case hoge =o 7B Teste -F (sc.o7 W/m °C) and Nu. eg PTO. 
Since T >. , the Nusselt number (and heat exchanger per- 
ave uniform 


formance) is lower when non-uniform flows exist. 
An enlarged view of the region of non-uniform flows is shown 


terres, 9. On Fig. 9 curve (c-d} at any AP level represents the 
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average flow rate with four tubes having a flow rate on curve (u-a) and 
one tube on curve (m1), the stable flow curves. Then SW ag = “Wie “ We 
Curves (s-f), (g-h) and (o-j) are similar curves for flow splits of 3-2, 
2-3 and 1-4 tubes on curves (u-a) and (m-1). As the total flow rate is 
increased, uniform flow exists through the peak at point(a) and on to 
point (b). At point (b) the pressure drop decreases suddenly to point 
(c) where one tube has a higher flow rate while the other four tubes 
have lower flow rates. The total flow rate is: 3Wy = WW = Aw ae 
As flow is further increased, flow increases in each of the tubes 
(along the pressure drop curve) with the average flow rate shown as 
line (c-d). Point (d) is the point depicted in Figs. 7 and 8. As flow 
is increased slightly above the flow corresponding to point (d), the 
pressure drop falls from points (d) to (e) and now two tubes have 
-higher flow rates and three tubes have lower flow rates. As the flow 
is increased, the flow in each tube increases with the average flow 
rate shown on the line (e-f). This process continues until uniform 
flow in all five tubes occurs at point (k). In the case of infinite 
tubes, the average flow rate would follow the line (a-1) with flow in 
the tubes corresponding to the flow at either point (a) - (lower flow 
rate) - (1) - (higher flow rate). The mix of tubes with higher and 
lower flow rates is determined by the average fluid flow rate. 
For the case of decreasing flow rate, non-uniform flows also 
occur. As flow is decreased to point (m) the pressure drop decreases 


with all tubes having uniform flow. The pressure drop suddenly in- 


creases from points (m) to (n) at which one tube has lower flow and 
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four tubes have higher flow. As flow is decreased further, flow in each 
tube decreases with the average flow rate given by line (n-o). Addi- 
tional tubes "pop-over" to lower flow at points (0), (qa), (s) and (c). 
All the tubes have the same flow rate at point (b). In the case of 
infinite tubes the average flow rate would follow line (m-u) with flow 
in the tubes corresponding to the flow at either point (u) - lower 

flow rate - or (m) - higher flow rate. 

Figure 10 shows the effects of non-uniform flow on the Nusselt 
number as a function of Reynolds number for the case of five tubes. 
The solid line is for increasing flow and the dashed line is for de- 
creasing flow. The Nusselt number is reduced significantly between 
points (b) and (k) for increasing flow and between points (m) and (b) 
for fecreasine flow. These are the regions of the curve where non- 
uniform flow exists. 

For the case of infinite tubes, the effect of non-uniform flow 
on exchanger performance is shown in Figs. 1l and 12. As the tube 
diameter is decreased, the region of non-uniform flow is extended 
to higher Reynolds numbers but the Nusselt number decreases less. 
Increasing the tube length has an effect on exchanger performance 


Similar to that of decreasing the tube diameter. 


b) Constant Head Source 
For the constant head source, as the pressure drop is in- 
creased beyond point (a) of Fig. 9, the average flow rate will in- 


crease to point (1) regardless of the number of tubes. With con- 





iS 


tinual increase in constant head, operation can not be sustained at 
flow rates between points (a) and (1). For the case of continual de- 
creasing constant head, all tubes remain operating along line (1) to 


(m). Then the:flow drops to point (u) in all tubes. 


c) Variable Head-Flow Pump 


- For a centrifugal pump source with the pump characteristic 
shown as lines (a-b) and (c-d) in Fig. 13, non-uniform flows can also 
exist. The situation is much the same as for the case of a constant 
head source except that the average flow follows the pump characteristic 
instead of the constant head line. As pump speed is increased above 
the point (a) all tubes shift to operation at point (b). The ex- 
changer can not operate between points (a) and (b) at the higher 


speed nor between points (c) and (d) at the lower speed. 


DISCUSSION 


Variations in the tube wall and inlet temperature can affect 
the values of exponents m and n (Figs. 2 and 3). Together with 
the u vs. T relation (Fig. 1), changes in exponents m and n 
can cause differences in the shape of the AP vs. w curve determined 
by Eq. (11). The AP a w curves for a single tube (L = 4’, 
D = 0.25") are plotted in Fig. 14 for different combinations of wall 


and inlet temperature. The tube with Ts = 40°F and es 300°F 


exhibits a region in which non-uniform flow can exist. If T, is 
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changed from 300°F to 200°F, there is only a slight dip in the curve 
and the region of non-uniform flow is narrower. If Ty is changed 
from 40°F to 100°F with T; = 300°F, no region of non-uniform flow 
exists. 

Results were similar for calculations done using another oil, 


Mobil automotive oil HD-30. 


CONCLUSIONS 


Laminar flow, parallel passage heat exchangers can experience 
maldistribution of flow and a reduction in the effective heat transfer 
coefficient in a range of Reynolds number which varies with length, 
tube diameter and the levels of the temperature - wall and inlet. 

For Mobil aviation oil 120 with the u vs. T relation of 
Fig. 1, the exponents m and n in the Baie correction ratios, 
Eqs. (3) and (4),are shown to vary with both wall and bulk temperatures, 


Figs. 2 and 3. 
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NOMENCLATURE 

s specific heat at constant pressure 

D tube diameter 

ia friction factor (variable viscosity) 

. friction factor (uniform viscosity evaluated at the bulk 
temperature) 

G mass flow rate 

hace effective heat transfer coefficient 

Ike thermal conductivity 

iG tube length 

m exponent for friction factor correction 

n exponent for Nusselt number correction 

ny number of tubes . 

Nu eg effective Nusselt number (variable viscosity) 

Nu Nusselt number (variable viscosity) 

Soe Nusselt number (uniform viscosity at the bulk temperature) 

Y pressure along the tube 

AP tube pressure drop 

¢ tube radius, radial direction 

Re, Reynolds number evaluated at Lh 

i temperature as function of r or x 

Bee average exit temperature in non-uniform flow condition 

T, bulk temperature 

T wall temperature 
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T, inlet temperature 

—. average exit temperature for uniform flow condition 
uniform 

u fluid velocity in x - direction 

W total flow rate 

W individual tube flow rate 

x direction along tube 


Greek Symbols 


0 fluid density 
u viscosity as function of T 
Hy viscosity evaluated at T, 


L viscosity evaluated at qT, 
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